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Dynamical correlations of various local operators are studied in the spin-half two-channel Kondo 
(2CK) model in the presence of channel anisotropy or external magnetic field. A conformal field 
theory-based scaling approach is used to predict the analytic properties of various spectral functions 
in the vicinity of the two-channel Kondo fixed point. These analytical results compare well with 
highly accurate density matrix numerical renormalization group results. The universal cross-over 
functions interpolating between channel-anisotropy or magnetic field-induced Fermi liquid regimes 
and the two-channel Kondo, non-Fermi liquid regimes are determined numerically. The boundaries 
of the real 2CK scaling regime are found to be rather restricted, and to depend both on the type 
of the perturbation and on the specific operator whose correlation function is studied. In a small 
magnetic field, a universal resonance is observed in the local fermion's spectral function. The 
dominant superconducting instability appears in the composite superconducting channel. 
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I. INTRODUCTION 

Deviations from Fermi liquid-like behavior observed 
e.g. in the metallic state of high-temperature cuprate 
superconductorSj-i^ or in heavy fermion systems'^''* 
prompted physicists to look for new non-Fermi liquid 
(NFL) compounds. So far a large number of such ex- 
otic compounds has been found and investigated. In 
these systems electrons remain incoherent down to very 
low temperatures and the usual Fermi liquid descrip- 
tion breaks down. To our current understanding, NFL 
physics may arise in many different ways: it can occur 
due to some local dynamical quantum fluctuations often 
described by quantum impurity modelSf^>^ii it can also 
be attributed to the presence of the quantum fluctua- 
tions of an order parameter or some collective modes, 
as is the case in the vicinity of many quantum phase 
transitions)^ or for the prototypical example of a Lut- 
tinger liquid)^ '^°i^^'^^i^^ where electrons are totally dis- 
integrated into collective excitations of the electron gas. 
NFL physics can also appear as a consequence of disorder 
like e.g. in disordered Kondo alloysj ^"^'^^ 

In this paper we study a variant of the overscreened 
multi-channel Kondo model: the spin-i, two-channel 
Kondo (2CK) model, which is the simplest prototypi- 
cal example of non- Fermi liquid quantum impurity mod- 
els. This model has first been introduced by Nozieres 
and Blandin^i^ and since has been proposed to de- 
scribe a variety of systems including dilute heavy fermion 
compoundSf§.tunneling impurities in disordered metals 
and doped semiconductors fiiii^ii^ More recently, the 
2CK state has been observed in a very controlled way 
in a double dot system originally proposed by Oreg and 
Goldhaber-Gordon3^ 

The two-channel Kondo model consists of a spin-i, 
local moment which is coupled through antiferromag- 
netic exchange interactions to two channels of conduc- 
tion electrons. Electrons in both channels try to screen 



the impurity spin. If the coupling of the spin to one of 
the channels is stronger than to the other then electrons 
in the more strongly coupled channel screen the spin, 
while the other channel becomes decoupled. However, 
for equal exchange couplings, the competition between 
the two channels leads to overscreening and results in a 
non- Fermi liquid behavior: Among others, it is character- 
ized by a non-trivial zero temperature residual entropy, 
a square root-like temperature dependence of the differ- 
ential conductance, a logarithmic divergence of the spin 
susceptibility and the linear specific heat coefficient at 
low temperaturesi^ This unusual and fragile ground state 
cannot be described within the framework of Nozieres' 
Fermi liquid theoryi^ 

Being a prototypical example of non-Fermi liquid 
models, the two-channel Kondo model (2CKM) has al- 
ready been investigated with a number of methods. 
These include non-perturbative techniques like the Bethe 
Ansatz, which gives full account of the thermodynamic 
properties r^i^ boundary conformal field theory;^ which 
describes the vicinity of the fixed points, and numeri- 
cal renormalization group (NRG) methodsi^ Further- 
more other less powerful approximate methods such as 
the Yuval-Anderson approach, ■^^ Abelian bosonizationj^ 
large-/ expansion^^^iM and non-crossing approximation^ 
have also been used to study the 2CKM successfully. 

Rather surprisingly, despite this extensive work, very 
little is known about dynamical correlation functions such 
as the spin susceptibility, local charge and superconduct- 
ing susceptibilities. Even the detailed properties of the 
T-matrix, essential to understand elastic and inelastic 
scattering in this non- Fermi liquid casCf^ have only been 
computed earlier using conformal field theory (which is 
rather limited in energy range) and by the non-crossing 
approximation (which is not well-controlled and is un- 
able to describe the Fermi liquid cross-over). '^^'■^^ It was 
also possible to compute some of the dynamical corre- 
lation functions in case of extreme spin anisotropy us- 
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ing Abelian bosonization results,-- though these calcu- 
lations reproduce only partly the generic features of the 
spin-isotropic model.— Local correlations in the Ander- 
son model around the non-Fcrmi liquid fixed point have 
already been investigated with the use of NRG, although 
in the absence of channel anisotropy and magnetic field4S 
However, a thorough and careful NRG analysis of the 
T = temperature T-matrix of the 2CKM has been car- 
ried out only very recentlyj^ii^ and the T 7^ analysis 
still needs to be done. 

The main purpose of this paper is to fill this gap by 
giving a comprehensive analysis of the local correlation 
functions at zero temperature using the numerical renor- 
malization group approach. However, in the vicinity of 
the rather delicate two-channel Kondo fixed point, the 
conventional NRG method fails and its further developed 
version, the density matrix NRG (DM-NRG)^ needs to 
be applied. Furthermore, a rather large number of mul- 
tiplets must be kept to achieve good accuracy. We have 
therefore implemented a modified version of the recently 
developed spectral sum conserving DM-NRG method, 
where we use non- Abelian symmetries in a flexible way 
to compute the real and the imaginary parts of various 
local correlation functions.— 

To identify the relevant perturbations around the NFL 
fixed point we apply the machinery of boundary confor- 
mal field theory. Then we systematically study how the 
vicinity of fixed points and the introduction of relevant 
perturbations such as a finite channel anisotropy or a fi- 
nite magnetic field influence the form of the dynamical 
response functions at zero temperature. We mainly fo- 
cus on the strong coupling regime of the 2CK model and 
the universal cross-over functions in the proximity of this 
region induced by an external magnetic field or channel 
anisotropy. We remark that these cross-over functions, 
describing the cross-over from the non-Fermi liquid fixed 
point to a Fermi liquid fixed point, as well as the re- 
sponse functions can currently be computed reliably at 
all energy scales only with NRG. However, we shall be 
able to use the results of boundary conformal field theory, 
more precisely, the knowledge of the operator content of 
the two-channel Kondo fixed point and the scaling di- 
mensions of the various perturbations around it, to make 
very general statements on the analytic properties of the 
various cross-over and spectral functions. 

We shall devote special attention to superconducting 
fluctuations. It has been proposed that unusual super- 
conducting states observed in some incoherent heavy 
fermion compounds could also emerge as a result of 
local superconducting correlations associated with two- 
channel Kondo physics i^i^i^ Here we investigate some 
possible superconducting order parameters consistent 
with the conformal field theoretical predictions, and find 
that the dominant instability emerges in the so-called 
composite superconducting channel^ as it was proposed 
by Coleman et al^ 

The paper is organized as follows. In Section [H] start- 
ing from the one-dimensional, continuum formulation of 



the 2CKM we connect it to a dimensionless approxima- 
tion of it suited to our DM-NRG calculations. We also 
provide the symmetry generators used in the conformal 
field theoretical and DM-NRG calculations. In Section 
mil we use boundary conformal field theory to classify the 
boundary highest-weight fields of the electron-hole sym- 
metrical 2CKM by their quantum numbers and identify 
the relevant perturbations around the 2CK fixed point. 
Based on this classification the fields are then expanded 
in leading order in terms of the operators of the free the- 
ory. In Section IIVI we describe the technical details of 
our DM-NRG calculations. In Sections IVl IVll and IVlIl we 
study the real and the imaginary parts of the retarded 
Green's functions of the local fermions, the impurity spin 
and the local superconducting order parameters. In each 
of these sections we first discuss the analytic forms of 
the susceptibilities in the asymptotic regions of the two- 
channel and single channel Kondo scaling regimes, as 
they follow from scaling arguments. Then we confirm 
our predictions by demonstrating how the expected cor- 
rections due to the relevant perturbations and the leading 
irrelevant operator present themselves in the DM-NRG 
data. Furthermore we determine the boundaries of the 
2CK scaling regimes and derive universal scaling curves 
connecting the FL and NFL fixed points for each oper- 
ator under study. Finally, our conclusions are drawn in 
Section 



II. HAMILTONIAN AND SYMMETRIES 

The two-channel Kondo model consists of an impurity 
with a magnetic moment S — ^ embedded into a Fermi 
liquid (FL) of two types of electrons (labeled by the flavor 
or channel indices a = 1,2), and interacting with them 
through a simple exchange interaction, 

+ EE V r ScUk)a,,c^,ik') . 

a ^ -I-Df J-Dp 

Here cj^, ^(k) creates an electron of flavor a in the / = 
angular momentum channel with spin /i and radial mo- 
mentum k measured from the Fermi momentum. In the 
Hamiltonian above we allowed for a channel anisotropy 
of the couplings, Ji 7^ J2, and denoted the Pauli ma- 
trices by tJ. In the flrst, kinetic term, we assumed a 
spherical Fermi surface and linearized the spectrum of 
the conduction electrons, ^(fc) « vpk = k, but these as- 
sumptions are not crucial: Apart from irrelevant terms in 
the Hamiltonian, our considerations below carry over to 
essentially any local density of states with electron-hole 
symmetry. The fields cjj^(fc) are normalized to satisfy 
the anticommutation relations 

[ci^k), c^^,(fc')} = 5^,p5^,u5 [k - k') , (2) 
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and therefore the couphngs Ja are just the dimensionless 
couphngs, usuaUy defined in the hterature. Since we are 
interested in the low-energy properties of the system, an 
energy cut-off D^^- is introduced for the kinetic and the 
interaction energies. In heavy fermion systems, this large 
energy scale is in the range of the Fermi energy, Dp ^ 
Ep, while for quantum dots, it is of the order of the single 
particle level spacing of the dot, Se or its charging energy, 
Ec, whichever is smaller. 

The Hamiltonian above possesses various symmetries. 
To see it, it is worth to introduce the left-moving fermion 
fields, 



dk e-''^^ 



(3) 



-Df 



and to rewrite the Hamiltonian as 
da; 



Then the total spin operators defined as 



da; 



(4) 

(5) 
(6) 



commute with the Hamiltonian and satisfy the standard 
SU(2) algebra. 



ijk 



(7) 



In the previous equations we suppressed spin indices and 
introduced the normal ordering : ... : with respect to the 
non-interacting Fermi sea. In a similar way we can define 
the "charge spin" density operators, for the channels a = 
1,2 as 

1 



- : Vl(a;)V'a(a;) 



C-{x) ^ VoT(^)V'ai(2;), C+{x) ^ ^i^(x)i^i^{x) , 
Ctix) ^ Cl{x)±iCl{x) , (8) 

and the corresponding symmetry generators 

dx 



27r 



(9) 



The generators which are related to the electron-hole 
symmetry;^ satisfy the same SU(2) algebra as the 



-s, 



(10) 



SUc2(2) X SU5(2) in the charge and spin sectors for ar- 
bitrary couplings, Ja- 

To perform NRG calculations, we use the following ap- 
proximation of the dimensionless Hamiltonian^^ 



Df{ 



+ E E ^»(/la,M/n+l,a.p + h.C.) , (11) 

n— a.^.y 

with A a discretization parameter and J a = 4Jc/(l + 
A~^). The operator, /o creates an electron right at the 
impurity site, and can be expressed as 



/o,( 



1 



Dp 



(12) 



The Hamiltonian Eq. (fTT|) is also called the Wilson chain: 
it describes electrons hopping along a semi-infinite chain 
with a hopping amplitude t„ ~ A^"/^, and interacting 
with the impurity only at site 0. In the NRG procedure, 
this Hamiltonian is diagonalized iteratively, and its spec- 
trum is used to compute the spectral functions of the 
various operators.— 

We remark that the Wilson Hamiltonian is not iden- 
tical to J^, since some terms are neglected along its 
derivatioufSE Nevertheless, similar to J^, the Wilson 
Hamiltonian also possesses the symmetry SU(7i(2) x 
SUc2(2) X SU5(2) for arbitrary Ji and J2 couphngs4^ 
The corresponding symmetry generators have been enu- 
merated in Table [J We can then use these symmetries to 
label every multiplet in the Hilbert space and every op- 
erator multiplet by the eigenvalues ^ ^ = j(j + 1) and 
'la'a — Ca{ca + !)• Throughout this paper, we shall use 
these quantum numbers to classify states and operators. 

In the presence of a magnetic field, i.e., when a term^ 



Hmagn = ~9^-bB 



(13) 



is added to the symmetry of the system breaks down 
to SUci(2) X SUc2(2) X U5(l), with the symmetry \Js{l) 
corresponding to the conservation of the z-component of 
the spin, (see Table H]). In the rest of the papers we 
shall use units where we set qijlb = 1- 



III. THE NON-FERMI LIQUID FIXED POINT 
AND ITS OPERATOR CONTENT 

For Ji ~ J2 ~ J and in the absence of an external mag- 
netic field, the Hamiltonian, possesses a dynamically 
generated energy scale, the so-called Kondo temperature. 



K 



Dp e 



and they also commute with the Hamiltonian, Eq. ([?]). 
Thus the Hamiltonian has a symmetry SUci(2) x 



The definition of Tk is somewhat arbitrary. In this paper, 
Tk shall be defined as the energy lu at which for Ji = J2 
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Symmetry group 




Generators 


oo 

SUc.(2) C+ = J2(^^)"fUrfL.i 




oo 

2 / J / J \ J n,(y.,fij n,a,iJ, ^ j 7 
n— f-i 


SUs(2) 


J 


oo 

- 9 + iV V a /■ 






n=Oa,f.L.h' 



TABLE I: Generators of the used symmetries for the two-channel Kondo model computations. Sites along the Wilson chain 
are labeled by n whereas a and fi, v are the channel and spin indices, respectively. 



the spectral function of the composite fermion drops to 
half of its value assumed at w = (for further details 
see the end of this Section and Fig. [T|). For B = and 
Ji = J2, below this energy scale the physics is governed 
by the so-called two-channel Kondo fixed point. 

The physics of the two-channel Kondo fixed point and 
its vicinity can be captured using conformal field the- 
ory. The two-channel Kondo finite size spectrum and its 
operator content has first been obtained using boundary 
conformal field theory by Affleck and Ludwig.— However, 
instead of charge SU(2) symmetries, Affleck and Ludwig 
used flavor SU(2) and charge U(l) symmetries to obtain 
the fixed point spectrum»2S The use of charge SU(2) sym- 
metries, however, has a clear advantage over the flavor 
symmetry when it comes to performing NRG calcula- 
tions: While the channel anisotropy violates the flavor 
symmetry, it does not violate the charge SU(2) symme- 
tries. Therefore, even in the channel anisotropic case, we 
have three commuting SU(2) symmetries. If we switch 
on a local magnetic field, only the spin SU(2) symmetry 
is reduced to its U(l) subgroup. Using charge symme- 
tries allows thus for much more precise calculations, and 
in fact, using them is absolutely necessary to obtain sat- 
isfactorily accurate spectral functions, especially in the 
presence of a magnetic field. 

To understand the fixed point spectrum and the oper- 
ator content of the 2CKM, let us outline the boundary 
conformal field theory in this SUci (2) x SUc2 (2) x SUs (2) 
language. First, we remark that the spin density oper- 
ators, J^{x) satisfy the 811(2)^=2 Kac-Moody algebra of 
level k = 2, 



[.rix),,Pix')] = -S^^S'ix~x') 

+ i2t: 5{x-x')e'^^J^{x) , (14) 

while the charge density operators, C^{x) defined in the 
previous section satisfy the Kac-Moody algebra of level 
fc = 1: 



Q(x),C^(x')J = -S'^S^^S'{x-x') 

+ i2TrSa,06{x-x')e'^''C^ix) 



struction to write the kinetic part of the Hamiltonian as 

J^fo = '^Cl + '^C2 + J^s + -^i , (15) 
-^Ca = 7; I ■ Ca{x)Ca{x) : , 



1 rdx 
37 2^ 



In Mq., the first two terms describe the charge sectors, and 
have central charge c = 1, while describes the spin 
sector, and has central charge c = 3/2. The last term 
corresponds to the coset space, and must have central 
charge c — 1/2, since the free fermion model has central 
charge c = 4, corresponding to the four combinations of 
spin and channel quantum numbers. This term can thus 
be identified as the Ising model, having primary fields 
l , fj, e with scaling dimensions 0, 1/16, 1/2, respectively. 
We can then carry out the conformal embedding in the 
usual way, by comparing the finite size spectrum of the 
free Hamiltonian with that of Eq. p5|) . and identifying 
the allowed primary fields in the product space. The 
the fusion rules obtained this way are listed on the left 
side of Table [TTl The finite size spectrum at the two- 
channel Kondo fixed point can be derived by fusing with 
the impurity spin (which couples to the spin sector only), 
following the operator product expansion of the Wess- 
Zumino-Novikov-Witten model, 1/2 (g) 1/2, 1/2 « 
1/2 ^ Oel, 1/2®1 1/2 (seeRHS of TableHH). Finally, 
the operator content of the fixed point can be found by 
performing a second fusion with the spin. The results of 
this double fusion are presented in Table [1111 In Table Hill 
the leading irrelevant operator, J71i</'s, is also included. 
Although it is not a primary fieldj^ close to the 2CK 
fixed point, this operator will also have impact on the 
form the correlation functions. 

What remains is to identify the scaling operators in 
terms of the operators of the non-interacting theory. In 
general, an operator of the non-interacting theory can 
be written as an infinite series in terms of the scaling 
operators and their descendants. Apart from the Ising 
sector, which is hard to identify, we can tell by looking 
at the various quantum numbers of an operator acting on 
the Wilson chain, which primary fields could be present 
in it. In this way, we can identify, e.g. (^a as the spin 
operator S. Thus the spin operator can be expressed as 



We can use these current densities and the coset con- 



s' = + 



(16) 
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TABLE 11: Left: Primary fields and the corresponding finite 
size energies at the free fermion fixed point for anti-periodic 
boundary conditions. States are classified according to the 
group SUci(2) X SUc2(2) x SUs(2) and the Ising model. The 
excitation energies Ehoc are given in units of 2-n /L, with L the 
size of the chiral fermion system. Right: Finite size spectrum 
at the two-channel Kondo fixed point. 



where the dots stand for all the less relevant operators 
that are present in the expansion of S, and some high- 
frequency portions which are not properly captured in 
the expansion above. The weight, Ag^ can be deter- 
mined from matching the decay of the spin-spin corre- 
lation function at short and long times. This way we end 
up with As ~ 1/ y/Tx- 

We remark that there are infinitely many operators 
that contain the scaling fields in their expansion. As 
an example, consider the operators (jy^^i- Here the label 

~ {Tii} refers to the spin components oi a, j = 1/2 
spinor, while t — ± refer to the charge spins of a charge 
c — 1/2 spinor. To identify the corresponding operator 
on the Wilson chain, we first note that /g x o- transforms 
as a spinor under spin rotations. It can easily be seen that 
the operator /q = i(Tj,/o,i also transforms as a spinor. 
We can then form a four-spinor out of these operators, 
'Ji = {fo I cry fo 1 a}- It is easy to show that 71 transforms 
as a spinor under SUci(2) rotations as well, thus (f>'^1 

could be identified as 71 = {/o,i,o-, /o.i^o-l- 
However, we can construct another operator, = 
/q ^S(7 and its counterpart, = iayFi, and form a four- 
spinor out of them: Ti = {^"i o-' ^1 o-i- This operator 
has the same quantum numbers as 71, and in fact, both 
operators' expansion contains 4'^^i- 

The operator cf)'^ is of special interest, since it is rel- 
evant at the two-channel Kondo fixed point, just like 
the spin. Its susceptibility therefore diverges logarith- 
mically. Good candidates for these operators would be 
X^crcr' ^<y<y'li'^l2 J since these are spin singlet operators 
that behave as charge 1/2 spinors in both channels. The 
T = t' = + component of this operator corresponds to 
the superconducting order parameter 



Osc — foA,1:fo,2,l foA,ifo,2,t ' 



(17) 



while the H — components describe simply a local opera- 
tor that hybridizes the channels, ^ fo 1 a-fo 2 a- 



Another candidate would be the operator, 
J2crcr' ^cr(T'ri'^72 ■ This Operator is also a local 
singlet, and has charge spins ci = C2 = 1/2. It 
contains the following component of the composite 
superconducting order parameter 



Oscc = fo.iSa i(Jyfl2 



(18) 



From their transformation properties it is not obvious, 
which one of the above superconducting order parame- 
ters gives the leading singularity. However, NRG gives 
a very solid answer and tells us that, while the suscepti- 
bility of the traditional operator does not diverge as the 
temperature or frequency goes to zero, that of the com- 
posite order parameter does. It is thus this latter oper- 
ator that can be identified as cj)'^ . Note that, in case 
of electron-hole symmetry, the composite hybridization 
operator 



(19) 



has the same singular susceptibility as Oscc since they 
are both components of the same tensor operator. This 
is, however, not true any more away from electron-hole 
symmetry. Furthermore, superconducting correlations 
are usually more dangerous, since in the Cooper chan- 
nel any small attraction would lead to ordering when a 
regular lattice model of two-channel Kondo impurities is 
considered. 

The knowledge of the operator content of the two- 
channel Kondo fixed point enables us to describe the 
effects of small magnetic fields and small channel 
anisotropics (Ji 7^ J2)- For energies and temperatures 
below Tk, the behavior of the model can be described 
by the slightly perturbed two-channel Kondo fixed point 
Hamiltonian. For Ji w J2 and in a small magnetic field, 
B <C Tk, this Hamiltonian can be expressed as 



n - rl2CK~ 

-F>Q Kq (pan 



DI'^ ho C^s+D^ XoJ^i 



-1/2 



(20) 



Here "H^ck 2CK fixed point Hamiltonian, and kq 

is the dimensionless coupling to the channel anisotropy 
field, 4>anis, whereas the effective magnetic field, ho, cou- 
ples to the "spin field" , 0s • Both of them are relevant 
perturbations at the two-channel Kondo fixed point and 
they must vanish to end up with the two-channel Kondo 
fixed point dX lo,T ^ 0. The third coupling, Aq, couples 
to the leading irrelevant operator (see Tab. IIIip . which 
dominates the physics when n = h = Q. The energy cut- 
off Do in Eq. ([^D]) is a somewhat arbitrary scale: it can 
be though of as the energy scale below which the two- 
channel Kondo physics emerges, i.e. Do ^ Tk- Then 
the dimensionless couplings kq, Aq and ho are approxi- 
mately related to the couplings of the original Hamilto- 



6 



Cl 


C2 


j 


/ 




scaling 

Upyi dLUl b 


correspondirif 


5 operators 








1 


1 


1 

2 








1 

2 





1 

2 




1 
2 




/ .4. 

71 = (/o,i,<T . 
Ti = (^0,1, <T 1 


(io-y^o,!)'^) 





1 

2 


1 

2 


a 


1 
2 


0^2 


72 





1 1 Q J 1 ^rr' I /d,l'^'^|'^!//d,2 " /(), 1 *5ct/o,2 

^ ^ ^ \ -/o,iO"y "S^CT crj,/J_2 -~fo,ii(ry Safo,2 

e i Cpania S {fl^^S fo ,x - fl^S fo^) 

1 I J-i^. ^(/d,ig/o,i +/d,2'?/o,2) 



TABLE III: Highest-weight operators and their dimensions x at the 2CK fixed point. Operators are classified by the 
symmetry group SUci(2) x SUc2(2) x SUs(2) and the scaling operators of the Ising model. The constants ci and C2 denote 
the charge spins in channels 1 and 2, respectively, while j refers to the spin, and I labels the scaling operators of the Ising 
model: fl , cr, e. Superscripts r, r' = ± refer to the two components of charge spinors, while a =|, J, label the components of a 
spin-±^ spinor. 



nian, Eq. ((TT|), as 

ho 
Ao 



B/Tk 
0(1) . 



J1-J2 

( Ji + 



(21) 

(22) 
(23) 



However, the arbitrary scale Dq in Eq. (I20p can be 
changed at the expense of changing the couphngs: Dq ~* 
D,Ko k{D), ho — > h{D) and Aq A(D) in such a 
way that the physics below Dq remains unchanged. This 
freedom translates to scaling equations, whose leading 
terms follow from the conformal field theory results, and 
read 



dK{D) 

dx 
dh{D) 

dx 
d\(D) 

dx 



1 



(24) 
(25) 
(26) 



with X = — log D. Solving these equations with the initial 
conditions, D — Dq ~ Tk and h = ho, k = kq, A = Aq, 
we can read out the energy scales at which the rescaled 
couplings become of the order of one, 



(Jl + J2)4 



K 



(27) 
(28) 



At these scales the couplings of the relevant operators 
are so large that they can no longer be treated as pertur- 
bations. Below T* the single channel Kondo behavior is 
recovered in the more strongly coupled channel, while 
can be interpreted as the scale where the impurity spin 
dynamics is frozen by the external field. 



The prefactors in Eqs. (|28p are somewhat arbitrary, 
and depend slightly on the precise definition one uses 
to extract these scales. In this paper, we shall use the 
spectral function of the composite fermion to define the 
scales Tk and T*. We define Tk to be the energy at 
which for Kn = the spectral function of the composite 
fermion takes half of its fixed point value (i.e. the value 
assumed at lo = 0). Whereas T* is the energy at which 
for Kji > it takes 75% of its fixed point value (see 
Fig. [ID. 

It is much harder to relate Th to a physically measur- 
able quantity. We defined it simply through the relation, 



Th = Ch—— 



Tk ' 



(29) 



where the constant was chosen to be Ch ~ 60. This way 
T/i corresponds roughly to the energy at which the NFL 
finite size spectrum crosses over to the low-frequency FL 
spectrum. 



IV. NRG CALCULATIONS 

Prior to discussing the analytic and numerical features 
of the response functions, let us devote this section to the 
short description of the NRG procedure used. All results 
presented in this paper refer to zero temperature. The 
NRG calculations were performed with a discretization 
parameter A = 2. The sum of the dimensionless cou- 
plings was J1 + J2 = 0.4 for each run. The NRG data were 
computed with a so-called flexible DM-NRG program^— 
which permits the use of an arbitrary number of Abelian 
and non- Abelian symmetries (see Tab. [l|, and incorpo- 
rates the spectral-sum conserving density matrix NRG 
(DM-NRG) algorithm.41 The DM-NRG method makes it 
possible to generate spectral functions that satisfy spec- 
tral sum rules with machine precision at T = tern- 
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FIG. 1: (color online) Spectral function Qp of the composite 
fermion operator, -Fo,i,T ^ ^ function of oj, and the definition 
of the scales Tk and T* . Tk is defined by the relation qf{uj — 
Tk,T = 0,Kr ^0) = ^Qf{lo = 0,T = 0,Kr ^ 0). For non- 
zero Kr the scale T* is defined through gri^^ = T* ,T — 
0,Kr) = ^gpiio ^ 0,T = 0,\Kr\). 



perature. For calculations with non-zero magnetic field 
the use of the DM-NRG method represents a great ad- 
vantage over conventional NRG methods,"*- which loose 
spectral weights and violate spectral sum rules. Conven- 
tional methods also lead to smaller or bigger jumps in 
the spectral functions at w = which hinder the compu- 
tation of the universal scaling functions provided by the 
scale The DM-NRG method solves all these prob- 

lems if a sufficient number of multiplets is kept. On an 
ordinary desk-top computer, however, we need to use as 
many symmetries as possible to keep the computation 
time within reasonable limits. 

In the present paper, where we study the electron- 
hole symmetrical case, it is possible to use the symmetry 
group SUci (2) X SUc2(2) x SUs(2) even in case of channel 
anisotropy. At these calculations the maximum number 
of kept multiplets was 750 in each iteration. This cor- 
responds to the diagonalization of « 85 matrices with 
matrix sizes ranging up to ~ 630, acting on the vector 
space of ~ 9000 multiplets consisting of « 106000 states. 
In the presence of magnetic field we used the symme- 
try group SUi:7i(2) x SUc2(2) x Us(l), and retained a 
maximum of 1350 multiplets in each iteration, that cor- 
responds to the diagonalization of w 150 matrices with 
matrix sizes ranging up to « 800 acting on the vector 
space of w 18000 multiplets consisting of « 73000 states. 

In the next sections, we shall see how the knowledge 
of the operator content of the two-channel Kondo fixed 
point can help us to understand the analytic structure of 
the various dynamical correlation functions obtained by 
NRG. 



V. LOCAL FERMIONS' SPECTRAL 
FUNCTIONS AND SUSCEPTIBILITIES 

Let us first analyse the Green's function of the local 
fermion, /q o- a 7« ■ composite fermion's (fJ g. q ^ 
^o,a) Green's function has already been looked into in 
detail in an earlier study of ours3i We shall therefore 
not discuss its analytic properties here but use it merely 
as a reference to define the various energy scales in the 
NRG calculations (see Fig.[T]). Let us note, however, that 
in the large bandwidth limit, u!,Tk «C Dp, the spectral 
function of the composite fermion and that of the local 
fermion are simply related. 



(30) 



Thus, apart from a trivial constant shift and a minus 
sign, the spectral function of the local fermion is that 
of the composite fermion, and all features of gp are also 
refiected in gf. 

Before we discuss the NRG results, let us examine what 
predictions we have for the retarded Green's function of 
the operator f^^.^ from conformal field theory. By look- 
ing at its quantum numbers, this operator can be identi- 
fied with the operator (pt^ (see Tab. IIIip. i.e. 



ft - 4 



f 



(31) 



with the prefactor Af oc 1/ ^Dp. Note that A/ is a com- 
plex number, it does not need to be real. The dots in the 
equation above indicate the series of other, less relevant 
operators and their descendants, which give subleading 
corrections to the correlation function of /, 



t 

0,fT,a ■ 



Further- 
more, the expansion above holds for the long time behav- 
ior. The "short time part" of the correlation function of 
/q Q, is not captured by Eq. , and gives a constant 
to Gf{uj) of the order of ^ l/Dp. Thus, apart from a 
prefactor Aj, a constant shift and subleading terms, the 

Green's function of fl^ ,^ is that of the field (j)'^- As 
we discuss it shortly in Appendix El the Fourier trans- 
form of the Green's function of any operator of dimension 
X — 1/2 is scale invariant around the two-channel Kondo 
fixed point. Since and thus /q ^ a have a scaling di- 
mension 1/2 at the 2CK fixed point, it follows that the 
dimensionless retarded Green's function, Dp Qf{ijj), is 
also scale invariant)^ 



Dp Gfiuo^T) 



D- 



dD 



= 9f 
= 0. 



LU T 

D'D 



,K{D),h{D),\{D),. 



(32) 



From Eq. p2p , we can deduce various important prop- 
erties. Let us first consider the simplest case, T = and 
K = h = 0. Then setting the scale D to Dq ~ Tx we 
have 



-K,h,T=0 

9f 



i^) = .9/ 



TT' "'■ 



(33) 
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Let us now rescale D — > \u)\, and use the fixed point 
scaling equation (pS)) to obtain A (I?), 



5/ 




(34) 



Assuming that this function is analytic in its second ar- 
gument we obtain for ^ Tk 



9f H 




(35) 



with g± f and g'_^ ^ some complex expansion coefhcients. 
Here the subscripts ± refer to the cases a; > and < 0, 
respectively. As we discussed above, the constants g± f 
depend also on the short time behavior of Qf{t), and are 
not universal in this sense. These constants are not inde- 
pendent of each other. They are related by the constraint 
that the Green's function must be analytic in the upper 
half-plane. Furthermore, electron-hole symmetry implies 
that g+f^g.f and g+ j = -{g'^ /)*• 

Relations similar to the ones above hold for the dimen- 
sionless spectral function. It is defined as 



1 



Qf{uj) = -- Img/(w) , 



(36) 




and assumes the following simpler form at small frequen- 
cies in case of electron-hole symmetry. 



(37) 



For Lo 3> Tk the scaling dimension of the local fermion 
is x^^'^'^ = 1/2 corresponding to an tj-independent spec- 
tral function. Perturbation theory in J amounts to log- 
arithmic corrections of the form: 1/2 — est/ log^(TR-/aj), 
as it is sketched in the upper parts of Fig.-s [2] and [31 

For T 7^ 0, and k, = h = Q using similar arguments as 
before, but now rescaling D T wc find 



9f i^) = 9f[j;,^Ao 




(38) 



Then by expanding gf we obtain the following scaling 
form for the low temperature behavior of the spectral 
function, 



0f (W) 




,(39) 



with Qf and Qf universal scaling functions. Note that 
we made no assumption on the ratio to/T, but both lu and 




log(t) logfe) 



r>o,A'B = o,s = o 




FIG. 2: (color online) (top) Sketch of the dimensionless spec- 
tral function gj = DpQf of /gj^, and (bottom) the real 
part of its dimensionless Green's function, Re gf = Df Re Gf 
for r > and Kr = 0,5 = as a function of log {uj/Tk)- 
Asymptotics indicated for u) < Tk were derived through scal- 
ing arguments. The large oj-behavior is a result of perturba- 
tion theory. 



T must be smaller than Tk- The asymptotic properties 
of 0/ and Q j can be extracted by making use of the 
facts that (z) gf{uj,T) must be analytic for lu -^T, (ii) 
that Eq. ([55)1 should reproduce the T ^ results in the 
limit uj ^ T, and {in) that by electron-hole symmetry, 
Qf must be an even function of uj. The issuing asymp- 
totic properties together with those of the other scaling 
functions defined later are summarized in Table [V] The 
asymptotic properties of the real part, Re gf, can be ex- 
tracted from those of by performing a Hilbert trans- 
form 



(40) 



with V the principal part. The obtained features are 
sketched in Fig. [2] for T > and k = /i = 0. 

Let us now investigate the effect of channel anisotropy, 
i.e. K ^ at T = temperature and no magnetic field 
h = 0. In this case, we can rescale D to D = \uj\ to obtain 



(w) 




(41) 
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FIG. 3: (color online) (top) Sketch of the dimensionless spec- 
tral function of /qi^,: = Dp Qf, and (bottom) the real 
part of its dimensionless Green's function: Re g/ — Dp Re Qf 
for T = 0,Kr > and B = as a function of log (lj/Tk)- 
Asymptotics indicated for u < Tk were derived through scal- 
ing arguments. The large oj-behavior is a result of perturba- 
tion theory. 



with T* the anisotropy scale defined earlier. The su- 
perscripts ± refer to the cases of positive or negative 
anisotropics: the superscript "+" is used when the cou- 
pling is larger in the channel where we measure the 
Green's function of /q ^ ^. The asymptotics of the univer- 
sal functions JC^ and K,^ can be obtained through similar 
scaling arguments as before and they differ only slightly 
from those of 0/ and 9/ (see Table fVl for a summary). 
The properties of ^J'''~'^(w) are summarized in Fig. [31 
A remarkable feature of the spectral function is that it 
contains a correction ^ yjT* /\ui\. This correction can 
be obtained by doing perturbation theory in the small 
parameter k[oj) at the two-channel Kondo fixed point. 

From the asymptotic forms in Table |V] we find that 
in the local fermion's susceptibility a new scale, TJ* ~ 

\/T*Tk appears as a result of the competition be- 
tween the leading irrelevant operator and the channel 
anisotropy)^ It is only in the regime TJ* < oj < Tk that 
the leading irrelevant operator determines the dominant 
scaling behavior of the local fermion's susceptibility, i.e., 
we observe the true two-channel Kondo physics. The 
expected properties of p/ and the real part of its dimen- 
sionless Green's function 5/ in the presence of channel 
asymmetry are summarized in Fig. [31 These analytic ex- 
pectations are indeed met by our NRG results. 



Fig. [41(a) depicts the spectral function of Jqi^ for 
several values of Kn as a function of uj/Tk on a loga- 
rithmic scale. The overall scaling is very similar to the 
one sketched in Fig. [31 except that the high tempera- 
ture plateau is missing; this is due to the relatively large 
value of Tk, which is only one decade smaller than the 
bandwidth cut-off. Figures [H (6 — e) are the numerical 
confirmations of the asymptotics stated. In all these fig- 
ures dashed straight lines are to demonstrate deviations 
from the expected behavior. In Fig. [11(6) we show the 
square root-like asymptotics in the 2CK scaling regime 
for the channel symmetric case. This behavior is a con- 
sequence of the dimension of the leading irrelevant op- 
erator as it has just been discussed. In Fig. [41(c) the 
same asymptotics is shown in the same region in case 

of a finite channel anisotropy, whereas below them Fig. 

r /2 

[31((i) demonstrates (l/i.^) -like behavior resulting from 
the relevant perturbation of the 2CK fixed point Hamil- 
tonian with channel anisotropy. In Fig. [H (e) the FL-like 
w^-behavior is recovered below T* , which is typical of 
fermionic operators in the ICK scaling regimes. 

In Figs. [51 (a — h) we show the universal scaling curves, 
/C* that connect the two-channel and single channel fixed 
points at low- frequencies as a function of lo/T*. They 
were computed from runs with negative and positive val- 
ues of Kji. This universal behavior is violated for values 
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FIG. 4: (color online) (a) Dimensionless spectral function of 
.fo,i,a- Qf{i-^) = Dp Qf(i-o) as a function of uj/Tk for different 
values of Kr. {b — e) Numerical confirmations of the low- 
frequency asymptotics derived through scaling arguments in 
Sec. [Vl Dashed straight lines are to demonstrate deviations 
from the expected y^-like (6— c), l/y^-like (d) and oj^-like 
(e) behavior. In plots (c - e) T* /Tk = 2.4 x 10"''. 
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Scaling Function 


Asymptotic Form 


Scaling Variable 


2CK 


a; < 1 , K X 


X 


Scaling Regime 



c^/T* Tf* < Lu , T}* oc VT-Tk 

uj/T T<Lo 



)6 iflO/ 2 aoO|aoo/|ll 1/2 



Qs (x) 0s° 2: , 0^ sgn(x) 

Qs (x) 0| X , sgn(x) lx|^/2 

W K^X, sgn{x) + " ^ r^,,^.. ,^.2^ ^l/3 

iCs (x) Kg sgn(x) Ixl^''^ , sgn(x) 



TABLE IV: Asymptotic behavior of the universal cross-over functions. At finite temperature, the boundary of the two-channel 
Kondo scaling regime is set by the temperature. At zero temperature, the various boundaries of the 2CK scaling regime derive 
from the competition between the leading irrelevant operator and the relevant perturbation. 
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FIG. 5: (color online) Universal collapse of the dimensionless 
spectral functions, gf — Dp Qf (with / in channel 1) to two 
scaling curves, K.^ as a function of ui/T* for positive (6) and 
negative (a) values of Ku. 



of higher than the highest ones shown in Fig.[5l where 
T* becomes comparable to T^. 

The real parts of the local fermion susceptibilities are 
plotted in Fig. [S] for several values of Kj^. They were ob- 
tained by performing the Hilbert transformations numer- 
ically. They should show a three-peak structure based on 
the analytic considerations (see Fig. [3]). There are two 
low-frequency peaks clearly visible, associated with the 
cross-overs at T* and Tr-. Furthermore there should be 
a non-universal peak at the cut-off. For relatively large 
channel anisotropics, where T* ~ T/f, the former two 
peaks cannot be clearly separated in Fig. [HI Also, due to 
the large value of Tr- ^ the band cut-off. Dp, the peak 
at w '--^ Tk and the smeared singularity at lo = Dp merge 



to a single non-universal feature in our NRG curves. 

Let us now turn to the effect of a finite magnetic field, 
B ^ ioT the case T — 0, Kb, = 0. As h and k scale 
the same way in the 2CK scaling regime, the argument 
concerning the k 7^ case can be repeated with minor 
modifications. Now, however, the spin SUs(2) symmetry 
is violated, and therefore the spectral functions oi 

and /o Q, I become different, and they are no longer even 
either. Nevertheless, due to particle-hole symmetry, they 
are still related through the relations 

PfA^^yT.K^h,...) = pf^i{uj,T,K,-h,...) . (42) 
We are thus free to choose the orientation of the magnetic 




FIG. 6: (color online) Real part of the dimensionless Green's 
function: Re gf = Dp Re Qf (with in channel 1) as a 
function oi uj/Tk for different values of Kr. From among the 
three peaks sketched in Fig. |3] only the two peaks around T* 
and Tk are shown. 
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FIG. 7: (color online) Top: Dimensionless spectral function 
of /o,i,T : Qf,t = Dp for different values of B as a func- 
tion oiuo/TK on linear scale. Bottom: Dimensionless spectral 
function of /o.i.t (a) and of /o,i,| ib) for different values of B 
as a function of u/Tk on logarithmic scale. 



field downwards. Then, after rescaling D \u}\ we get 




(43) 



where the label a refers to the different spin components 
and Bf^a and Bj^a are yet another pair of universal cross- 
over functions. The asymptotic properties of the func- 
tions Bf^a and Bf\a are summarized in Table [V] 

Fig. [7] shows the spectral functions gf^a as a function 
of lo/Tk on linear and logarithmic scales for different 
magnetic field values. The same curves are depicted as 
a function of lo /Th in Fig. [SI which demonstrates the ex- 
istence of the universal scaling curves, Bf^^, i.e. that by 
using the scale, Th the local fermion's spectral functions 
can be scaled on top of each other for small enough mag- 
netic fields. In this magnetic field region, we find a peak 
at Th for the spin-f component of while at the same 
place there is a dip for the spin-| component. This is a 
remarkable feature that is associated with inelastic scat- 
tering off the slightly polarized impurity spin. In fact. 



3, 



0.25 




0.25 



FIG. 8: (color online) Universal collapse of the dimensionless 
spectral functions: ^/,f = Dp Qf^i and Qf^i = Dp Qf,i to two 
scaling curves: and Bf,i for sufflciently small, non-zero 

values of B as a function of uo/Th- 



the same uinversal features also appear in the spectral 
functions of the composite fermions, which we compute 
independently and which are directly related to those of 
the conduction electrons by Eqn. ([30]) .— The rescaled 
spectral functions QF,a{^) are shown in Fig[51 

Although this numerical evidence can be obtained by 
conventional NRG methods not using the density matrix, 
this is no longer true for the sum of the local fermions 
spectral function over the different spin components. In 
fact, for this quantity universal scaling curves in the pres- 
ence of magnetic field cannot be obtained using NRG be- 
cause of the increase in the size of the numerical errors at 
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FIG. 9: (color online) Universal collapse of the dimensionless 
spectral functions of the composite fermion operator: gp^i = 
Dp Qf^-^ and gp^i — Dp gp^i to two scaling curves: Sf,t and 
Bp^ l for sufficiently small, non-zero values of i? as a function 
of'uj/Th. 
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FIG. 10: (color online) Sum of the dimensionless spectral 
functions: ^/,| — Dp Qf,\ and = Dp gf.^ for different 
values of -B as a function of uj/Tk- 



low-frequencies and the mismatch between the positive 
and negative frequency parts of the spectral functions. 
The sum of the local fermion's spectral function over the 
two spin components is depicted in Fig.[Tn]as a function 
of uj/Tk- Here the splitting of the Kondo resonance in 
the energy-dependent scattering cross section appears as 
a minimum at a; ~ T^. Unfortunately, for even smaller 
magnetic fields the accuracy of our numerical data is in- 
sufficient to tell if the splitting of the Kondo resonance 
survives in the limit _B ^ 0, as conjectured in Ref. (23. In 
the data with B/Tk > 1.1 x 10~^, there seems to be al- 
ways a shallow minimum in the spectral function, and we 
see no indication for crossing of the curves as the magni- 
tude of the field is reduced. If there is indeed no crossing 
of the spectral functions and if the deviation from the 
•\/|w|-behavior indeed starts at tu ~ Th ^ B'^/Tk, which 
is the only natural assumption, then, from exact Bethe 
Ansatz results it would immediately follow that there 
must always be a splitting of the Kondo resonance, since 

EAi^fA^ = 0,5) - QfA^ - 0,0)] ^ bhtk/b)P. 

while EAqM^ = Th,B) - gf„{uj = 0,0)] - \B\ would 
follow from the pure -y/ 1 w | -dependence of the spectral 
function at B = 0. However, these analytical arguments 
do not constitute a real proof. 

With small modifications, the analysis presented in 
this subsection carries over to essentially any fermionic 
operator that has quantum numbers ci = j = 1/2 or 
C2 = j = 1/2 and has a finite overlap with the pri- 
mary fields (p^i and 0^2, only the high-frequency be- 
havior {lo > Tk) and the normalization factors become 
different. Typically, a local operator having the same 
charge and spin quantum numbers as (/)^q, will have a 
finite overlap with them. However, in some cases the in- 
ternal Ising quantum number of an operator may prevent 
an overlap and, of course, one can also construct opera- 
tors by, say, differentiating with respect to the time, that 
would correspond to descendant fields. 



VI. SPIN SPECTRAL FUNCTIONS AND 
SUSCEPTIBILITIES 

In this section, we shall discuss the properties of the 
spin operator, S, which is the most obvious example of a 
bosonic operator of spin j = 1 and charge quantum num- 
bers ci = C2 = that overlaps with the scaling operator 
(ps- There are, however, many operators that have the 
same quantum numbers: two examples are the so-called 
channel spin operator. 



Sc 



or a composite channel spin operator 

Sec = fI^^SJo^i - -Fo,2'?/o,2 



(44) 



(45) 



Our discussion can be easily generalized to these opera- 
tors with minor modifications. 

The analysis of the spin spectral function goes along 
the lines of the previous subsection. First we recall that 
the field 0s appears in the expansion of the spin operator. 



5 = As 



(46) 



with As ^ \/\JTk ^ l/V^o- Therefore, the appropriate 
dimensionless scale invariant Green's function (usually 
referred to as the dynamical spin susceptibility) is defined 
as 



9s \ j^,j^,Ko,ho,. 



= Tk Gs{oj,T,ko,...,Do) . 



(47) 



We shall not repeat here all the steps of the deriva- 
tion, only summarize the main results. In the absence 
of a magnetic field, h = the spectral function of the 
spin operator is odd. Furthermore, at T = and for no 
anisotropy, n — 0, the spectral function has a jump at 



- n 5 



0^ 



0S (^) 



sgn(w) 




(48) 



This jump corresponds to a logarithmically divergent 
dynamical susceptibility, Kexsi'^) = —^gQs{u>) cx 
HTK/io)/TK. 

For Lu 3> Tk the impurity spin becomes asymptotically 
free, decoupled from the conduction electrons, therefore 
its cij-dependence is set by its scaling dimension at the 
free fermion fixed point where Xg^'^'' = 0. It has the im- 
plication that its correlation fuction decays as a;~^ cor- 
responding to the Curie- Weiss susceptibility with loga- 
rithmic corrections present, known from Bethe Ansatz 
results and from perturbation theory. 

At finite temperatures T 0, but for k ~ h ~ 0, we 
obtain the following scaling form for T,uj Tk- 



9S (^) = '^S[7^) 




(49) 
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FIG. 11: (color online) Top: Sketch of the dimensionless spec- 
tral function of S: gs = TkQs = -Txlm xs(a;)/7r for T > 
and Kr = 0, B — as & function of \og{uj/TK)- Bottom: 
Sketch of gs = TkQs = -Txlm xs{uj) /-ir for T = and 
Ku 0, B — as a, function of log (uo/Tk)- Asymptotics in- 
dicated for uj < Tk were derived through scaling arguments. 
The large tj-behavior is a result of perturbation theory.— 

The asymptotic properties of the scaling functions Qg 
and 8s are listed in Table IVl 

In case of finite channel anisotropy but zero tempera- 
ture we obtain for uj <C Tk the scaling form 

fr-\.)^Ks[^y^^ts[^y.... (50) 

The asymptotic properties of ICs,)Cs a-re only slightly dif- 
ferent from those of 65, Qs (see Table IV)) : below T* the 
spectral function displays analytic behavior, while the 
regime w > T* is governed by non-analytical corrections 
associated with the 2CK fixed point. In this regime a fea- 
ture worth mentioning is the appearance of a correction, 
~ T*/lj to K-s, more precisely, the lack of a ^/\T*/uj\ 
correction. This is due to the fact that the anisotropy 
operator is odd, while the spin operator is even with re- 
spect to swapping the channel labels. Therefore there is 
no first order correction to the spin-spin correlation func- 
tion in K, and the leading corrections are only of second 
order, i.e., of the form k'^/uj. From the comparison of 
the terms in /Cg and ICs it also follows the existence of 




FIG. 12: (color online) Sketch of the real part of the di- 
mensionless Green's function of S, Re gs ~ Tk Re Xs{^) = 
Tk Re Gsico) for T, T* > as a function of log (lo/Tk). 

another cross-over scale, 

T** {T*^Tk\ , (51) 

that separates the regimes governed by the leading rel- 
evant and leading irrelevant operators. Here we used 
the subscript s to indicate that this scale T** is differ- 
ent from the scale Tj* introduced in relation to the local 
fermion's spectral function. The asymptotic properties 
of Qs oc xs{^^) for T>Q,Kr = and T = 0, Kr ^ are 
sketched in the upper and lower parts of Fig. 1111 while 
the behavior of the real part is presented in Fig. [T2l 

The expectations above are indeed nicely born out by 
the NRG calculations: Fig. [13] shows the impurity spin 
spectral functions as a function of w/Tr- for various Kr- 
s and their asymptotic properties. First, in Fig. 1131 (6) 
we show a very small logarithmic w-dependence that we 
observed below Tk at the 2CK fixed point. The am- 
plitude of this log(w)-dependence was reduced as we in- 
creased the number of multiplets. It appears that this 
behavior is not derived from the lognormal smoothing 
of the NRG data, and it may be due to some approxi- 
mations used in the spectral sum-conserving DM-NRG 
procedure. In Fig. [131 (c) we show the square root-like 
behavior around the 2CK Kondo fixed point which is 
attributed to the leading irrelevant operator, while Fig. 
[131 (d) shows that first order corrections coming from the 
scaling of the channel anisotropy are indeed absent just 
as we stated above, and only second order terms appear, 
resulting in an l/w-like behavior. Finally, Fig. [T3l(e) 
demonstrates the linear w-dependence, which is charac- 
teristic of most bosonic operators in the proximity of an 
FL fixed point. All these findings support very nicely the 
analytical properties summarized in Table [VI 

The spin spectral functions also collapse to a univer- 
sal scaling curve describing the cross-over from the two- 
channel Kondo to the single channel Kondo fixed points, 
when they are plotted against lo/T** . This universal data 
collapse is demonstrated in Fig. [U where the impurity 
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FIG. 13: (color online) (a) Dimensionless spectral function 
of S: gs ~ TkQs = — TkIhi xs(i^)/7r as a function of lo/Tk 
for different values of Kji- (b) Minute log (Lj)-dependence at 
the lowest frequencies diminishing as a function of the num- 
ber of kept multiplets. (c — e) Numerical confirmations of the 
low-frequency asymptotics derived from scaling arguments in 
Section IVII Straight dashed lines are to demonstrate devia- 
tions from the expected v^-like (c), w-like (d), and 1/ co-like 
behavior (e). In plots (d - e) T* /Tk = 7 x 10"^. 



spin spectral functions are plotted for various Kji val- 
ues. The data collapse works up to somewhat higher 
anisotropy values than for the local fermions' spectral 
functions as it is indicated by the i^/j-dependence of the 
scales r;* and Tf. 

The real part of the spin susceptibility was obtained 
through numerical Hilbert transformation, and is shown 
in Fig. [15] as a function of lo/Tk for various values of 
Kn. These curves meet the expected behavior sketched 
in Fig. [T2I they display a logarithmic increase at high- 
frequencies and saturate at values that correspond to 
Re xs ~ HTk/T*) I Tk. 

Let us finally discuss the case, T ~ n — Q but h ^ 
0. Then the components of S are distinguished by the 
magnetic field: The spectral function of S^- has almost 
the same features as for finite channel anisotropics. Since 
is a hermitian operator, its spectral function remains 
odd and acquires the following corrections in the different 
scaling regimes 



with the scaling functions Bs,z,l3s,z having the asymp- 
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FIG. 14: (color online) Universal collapse of the dimensionless 
spectral function of S: gs = Tkqs to the scaling curve, Ks 
as a function of uj/T* for sufficiently small, non-zero values of 
Kr. 



totic properties listed in Table IVl 

Note that in this case the first order correction coming 
from the magnetic field does not vanish, and leads to the 
appearance of a cross-over scale ~ y/T^TK- 

The perpendicular components of the impurity spin 
have somewhat different properties. First of all, the op- 
erators are not Hermitian, and therefore their spectral 
functions are not symmetrical. The spectral functions of 
the operators 5^ and are, however, symmetrical, and 
their Green's functions (and susceptibilities) are related 
through 
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(53) 



The corresponding dimensionless spectral functions, 
and gf as computed by our DM-NRG calculations are 
shown in Fig. [12] as a function of lo/Tk, while the uni- 
versal scaling with Lo/Tfi is confirmed for low-frequencies 
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FIG. 15: (color online) Real part of the dimensionless Green's 
function (susceptibility) of S, Re gs = —Tk Re xs{lo), as a 
function of uj/Tk, for different values of Kr. 
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FIG. 16: (color online) (a) Dimensionless spectral function 
of S+: Qs,+ = TkQs,+ , (b) of S^: gs,. = Tkqs,. and (c) of 
S~: Qs,- ~ TkQs,- for different values of _B as a function of 
uj/Tk- 



FIG. 17: (color online) Universal collapse of Qs^+ = TkQs,+ , 
Qs.z = TkQs,z and ^s,- = TkQs,- to the three scaling curves: 
Ss,+ ,Ss,2 and Bs.~ for sufficiently small, non-zero values of 
B as a function of uj /Th ■ 



in Fig. [171 This scaling also turned out to be valid for 
values of B higher than the ones for fermions (see Fig. 
[T7)) . The scaling functions Bs.z and Bs,± behave very 
similarly. This is somewhat surprising, since the naive 
expectation would be to have a resonance in Bs,+, just 
as in the local fermion's spectral function, that would 
correspond to a spin-flip excitation at the renormalized 
spin splitting, Th- 
excitations" . 

However, quite remarkably, a resonance seems to ap- 
pear in Xs zi'^)/^ ^ frequency uj ^ Th, while wc find 
no resonance in Xs ±{'^)/^- This can be seen in Fig. [T8l 
where T'^g{uj)/uj is plotted for the different spin compo- 
nents as a function of oj/Tk for various magnetic field 
values. This seems to indicate that the spin coherently 
oscillates between the spin up and spin down compo- 
nents, while its x,y components simply relax to their 
equilibrium value. 



VII. SUPERCONDUCTING CORRELATIONS 

In the last section, let us investigate the local super- 
conducting correlation functions. These deserve special 
attention, since many heavy fermion compounds display 
exotic superconducting phases that may possibly be in- 
duced by local two-channel Kondo physics The most 
obvious candidates for the corresponding local operators 
have been identified in Section llllj and are the local 
channel-asymmetric superconducting operator, Osc = 
/d 1 T'^o 2 X ~ /d 1 x/d 2 t' ^'^'^ composite fermion su- 
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FIG. 18: (color online) (a) q{ijj)/u: of 5*+: TkQs,+/'2u = 
tIqs,+ /2lo, (b) of S*": TkOs^J^ = TI,qs,./lo and (c) of S": 
TkQs -/2tt; = T^gs.- /2ijj for different values of i? as a func- 
tion of u/Tk. 

perconductor field, Oscc = /o I'S'f? 'io-yfo2- 

For the composite superconductor we find the expan- 
sion, 

Oscc = Ascc + • • • (54) 

where the expansion coefficient Ascc can be estimated 
from the high-frequency behavior of the correlation func- 
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tion up to logarithmic prefactors as Ascc ~ VTk/Dp. 
While for the impurity spin, one can exclude logarith- 
mic corrections to the expansion coefficient As in Eq. [351 
based upon the exact Bethe Ansatz results, this is not 
possible for the superconducting correlation function. In 
fact, we know that in the expansion of the composite 
fcrmion itself the correct prefactor is ^ J/VTk ''^ 
1/ {\/Tk \n{Dp/TK)) ■'^^ Therefore, similar logarithmic 
factors could appear in the prefactor Ascc- Neverthe- 
less, in the following, we shall disregard possible logarith- 
mic corrections, and define the normalized dimensionless 
and scale-invariant correlation function through the rela- 
tion. 



D2 



T, 



K 



Qscc{^^) = 5scc(w) 



(55) 



Apart from its overall amplitude and its high-frequency 
behavior, in the low-frequency scaling regimes the spec- 
tral function of the composite superconductor operator 
behaves the same way as that of (see Tab. lIIip . There- 
fore we merely state its asymptotics without further ex- 
planation. 

In the absence of anisotropy and magnetic field, k = 
/i = 0, for oj <^ Tk the spectral function becomes a uni- 
versal function, pscci'-^/T), whose behavior is described 
by the scaling form. 



Oscc (^) ~ Q 



sec 



^ ©^C7C (^) 



(56) 



while in the presence of anisotropy, but at T = tem- 
perature and for h = Q, the spectral functions behave 
as 



9scc (^) 



sec 




sec 



(57) 



Finally, in a finite magnetic field but for k — 
anisotropy and T — temperature the spectral function 
assumes the following scaling form. 



l^'scc 



Bscc I — 




(58) 

The properties of the the various scaling functions de- 
fined above are identical to those of the corresponding 
spectral functions of the , which were detailed in Ta- 
ble |Vl therefore they have not been included in Table |Vl 
The asymptotic properties are nicely confirmed by our 
NRG calculations. The dependence on the anisotropy, 
together with the ~ \/M , the ^ 1/uj and the ^ oj scaling 
regimes are plotted in Fig. [191 Here the high-frequency 
region, lo > Tr-, is also displayed, where the spectral 
function is roughly linear in the frequency, as dictated 
by the free fermion fixed point. 
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FIG. 19: (color online) (a) Dimensionless spectral function, 
Qscc ~ D%/Tk Qscc of the operator Oscc, as a function of 
oj/Tk for different values of Kr. (6). The very weak log [uj)- 
dependence at the lowest frequencies. This dependence is sup- 
pressed as we increased the number of kept multiplets. (c — e) 
Numerical confirmations of the low- frequency asymptotics de- 
rived from scaling arguments in Section rVIII Dashed straight 
lines are to demonstrate deviations from the expected ^/uo- 
like (c), cij-like (d) and l/ti;-like (e) behavior. In plots {d— e) 
T*/Tk = 7 X 10"^ 



Fig. [201 displays the real part of the dimensionless 
Green's function, that is essentially the real part of the 
superconducting susceptibility. This diverges logarithmi- 
cally for T* = 0, but for finite r*'s it saturates, corre- 
sponding to a susceptibility value 



Xscc 



Tk , (Tk 



Notice that there is a small prefactor in front of the log- 
arithm that arises from the asymptotically free behavior 
at large frequencies. 

The universal collapse of the low-frequency part of the 
curves in terms of cu/T* is shown in Fig. [^H The cross- 
over curve, JCscc (^) is very similar to the spin cross- 
over function, K-s, and displays a plateau at large fre- 
quencies from which it deviates as 1/w, until it finally 
reaches the linear frequency regime below T* . 

Application of a magnetic field has effects very similar 
to the anisotropy, as shown in the upper part of Fig. [HI 
In Fig. [22l the small logarithmic increase at small fre- 
quencies is more visible. As mentioned before, this in- 
crease is most likely an artifact of the spectral sum con- 
serving approximation of Ref. |4ll and it is due to the 
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FIG. 20: (color online) Real part of the dimensionless Green's 
function of Oscc- Re gscc ~ DI/Tk Re Gscc as a function 
of u)/Tk for different values of Kr. 



way this method redistributes spectral weights. This is 
based on the observation that the slope of the logarithm 
gets smaller if we increase the number of multiplets kept. 
These curves also collapse to a single universal curve as 
a function of a;/T^, as shown in the lower part of Fig.[22l 

Finally, in Fig. I23i we show the numerically ob- 
tained spectral function and the corresponding dimen- 
sionless susceptibility of the non-composite superconduc- 
tor, Osc = /d,i,T/o,2,i~/d,i,i'^d,2,T- Clearly, this spectral 
function displays no plateau below T/f , but it exhibits a 
linear in lj behavior below Tk, and correspondingly, the 
susceptibility Re xsc remains finite for cj ^ even in 
the absence of anisotropy and an external magnetic field, 
i.e., at the 2CK fixed point. 

This implies that, although its charge and spin quan- 
tum numbers would allow it, the expansion of this op- 
erator does not contain the scaling operator (j)^^ . This 
may be due to the difference in the Ising quantum num- 
bers, which we did not identify. Thus the dimension of 
the highest-weight scaling operator that appears in the 




FIG. 21: (color online) Universal collapse of gscc to the scal- 
ing curve tCscc as a function of ujjT* for sufficiently small, 
non-zero values of Kr. 
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FIG. 22: (color online) Top: Dimensionless spectral func- 
tion Qscc = -{D%I-kTk) Im xscc = -{D^I-rTk) Im ^scc 
of the composite superconductor operator Oscc for different 
values of B, as a function of uj /Tk- Bottom: Universal col- 
lapse of Qscc to the scaling curve Bscc for sufficiently small, 
non-zero values of B as a function of to /Th . 



expansion of Osc is a; = 1 and not 1/2, as one would 
naively expect based upon a simple comparison of quan- 
tum numbers. Turning on a small anisotropy or magnetic 
field does not influence substantially the spectral proper- 
ties of the corresponding Green's function, either. 



VIII. CONCLUSIONS 

In the present paper we gave a detailed discussion of 
the spectral properties of the two-channel Kondo model. 
We analyzed the properties of the correlation functions 
of various local operators in the presence of a channel 
anisotropy and an external magnetic field. In particular, 
we studied numerically and analytically the correlation 
functions of the local fermions, /a.o- = fo,a.a, the compo- 
nents of the impurity spin, S, the local superconductivity 
operator, Osc = fii'^yfl,, and the composite supercon- 
ductor operator, j\Sa iuyf\. The selection of these op- 
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FIG. 23: (color online) (a) Dimensionless spectral function 
of Osc- Qsc = Df qsc as a function of lu/Tk for different 
values of Kr, and (b) the real part of its dimensionless Green's 
function: Re gsc ~ Df Re Qsc- 




erators was partially motivated by conformal field theory, 
which tells us the quantum numbers and scaling dimen- 
sions of the various scaling operators at the two-channel 
Kondo fixed point There are, however, many operators 
that have quantum numbers identical with the scaling 
fields. Here we picked operators having the right quan- 
tum numbers, and at the same time having the largest 
possible scaling dimension at the free fermion fixed point, 
where Ji, J2 0. These are the operators, whose spec- 
tral functions are expected to have the largest spectral 
weight at small temperatures (among those having the 
same quantum numbers), and which are therefore the 
primary candidates for an order parameter, when a lat- 
tice of 2CK impurities is formed, as is the case in some 
Uranium and Cerium-based compounds. The operators 
above are, of course, also of physical interest on their 
own: the spectral function of fa^a is related to the tun- 
neling spectrum into the conduction electron see at the 
impurity site, the Green's function of S is just the dy- 
namical spin susceptibility that can be measured under 
inelastic neutron scattering, and finally the local super- 
conducting operators are candidates for superconducting 
ordering in heavy fermion materials. We remark that, in 
the electron-hole symmetrical case, the other components 
of the operator multiplet that contains the composite su- 
perconducting order parameter Oscc would correspond 
to a composite channel-mixing charge density ordering. 
Of course, the susceptibilities of this operator has the 
same properties as that of xs'cc(w)- 

In addition to these operators, there are two more 
operators of possible interest: the so-called composite 
Fermion's Green's function is related to the T-matrix, 
T{uo) that describes the scattering properties off a two- 
channel impurity (or the conductance through it in case 
of a quantum dot), and was already studied in detail in 
Refi^. A further candidate is the channel anisotropy op- 
erator. This has also a logarithmically divergent suscep- 
tibility, and would also be associated with a composite 
orbital ordering in case of a two-channel Kondo lattice 
system. However, the spectral properties of this latter 
operator are so similar to those of the composite super- 
conductor that we have decided no to show data about 
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FIG. 24: (color online) Top: Sketch of the various 2CK scal- 
ing regimes in the presence of channel anisotropy for the lo- 
cal fermions bounded by T" from below and for the spin 
bounded by T** , the crossover scale T** is also indicated. 
Bottom: Sketch of the 2CK scaling regime for the suscepti- 
bilities of the highest- weight fields bounded by T" , and the 
crossover scale Th besides. 



them. 

For the numerical calculations we used a flexible DM- 
NRG method, where we exploited the hidden charge 
SU(2) symmetries'*^ as well as the invariance under spin 
rotations to obtain high precision data. To identify 
the scaling operators in this case, we reconstructed the 
boundary conformal field theory of Affleck and Ludwig 
for this symmetry classification. We then established 
the scaling properties of the various dynamical correla- 
tion functions and identified the corresponding univer- 
sal cross-over functions and their asymptotic properties, 
based upon simple but robust scaling arguments. In this 
way, universal scaling functions describing the cross-over 
from the two-channel Kondo fixed point to the single 
channel Kondo fixed point (for Ji 7^ J2) and to the mag- 
netically polarized fixed point (for B ^ 0) have been 
introduced, which we then determined numerically. We 
emphasize again that presently these universal cross-over 
functions can only be determined through the applica- 
tion of DM-NRG, and in fact, for the scaling curves in 
the presence of a magnetic field the application of the 
DM-NRG method was absolutely necessary. 

Our numerical calculations confirmed all our analyt- 
ical expectations, and they confirmed that actually, in 
the presence of an applied magnetic field, or channel 
anisotropy, the two-channel Kondo scaling regime is 
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rather restricted, and it may also depend on the physical 
quantity considered. In FiglM]we sketched the regimes 
where the pure two-channel Kondo behavior can be ob- 
served. Notice that in the presence of anisotropy the 
two-channel Kondo scaling regime of the spin suscepti- 
bility has a boundary that differs from the boundary of 
the two-channel Kondo scaling regime of the T-matrix. 

Some of the spectral functions show rather remarkable 
features: In a magnetic field, e.g., the spectral function 
of the composite fermion, ^ shows a universal peak 
at a frequency lo = T^. This peak corresponds to spin- 
flip excitations of the impurity spin at the renormalized 
magnetic field. Remarkably, this peak is accompanied 
by a dip of the same size at the same frequency for spin 
down electrons. This dip is actually very surprising and 
is much harder to explain. Similar features appear but 
with opposite sign in the local fermions' spectral func- 
tions. Even more surprisingly, this resonant feature is 
completely absent in the spectral function of the spin op- 
erators, S^. 

One of the interesting results of our numerical anal- 
ysis was that only the composite superconductor Osc 
has a logarithmically divergent susceptibility. This is 
thus the primary candidate for superconducting ordering 
for a 2CK lattice system. We remark here that while 
for a single impurity the superconducting susceptibil- 
ity seems to have a rather small amplitude. Re xsc ~ 
Tk / Dpln{TK /T), in a lattice model the mass of the 
carriers is also renormalized, and therefore the band- 
width is expected to get renormalized as Dp As 
a result, the corresponding susceptibility can be rather 
large, and drive, in principle, a superconducting instabil- 
ity. Interestingly, although the results are still somewhat 
controversial)^ in the two-channel Kondo lattice these 
local superconducting correlations do not seem to induce 
a superconducting transitional This may be, however, an 
artifact of the standard two-channel Kondo lattice model, 
which does not account properly for the orbital and band 
structure of an /-electron material4^ We believe, that in 
a more realistic lattice of two-channel Kondo impurities a 
composite superconducting order develops, similar to the 
one suggested in Ref. S. However, DMFT + DM-NRG 
calculations would be needed to confirm this belief. 
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APPENDIX A: SCALING PROPERTIES OF 
TWO-POINT FUNCTIONS 

In this appendix, we discuss the scaling properties of 
various scaling functions. Essentially, we use the gener- 
alized Callan-Symanzik equations. For the sake of sim- 
plicity, let us first focus on the retarded Green's function 
of the z-component of the operator (j)s , 



-^{Ws[t)Alm )n e{t) 



(Al) 



and its Fourier transform, Q{u}, T). Let us investigate the 
scaling properties of this function in the absence of mag- 
netic field. From the fact that (ps is the field conjugate 
to the external "magnetic field" , h, and that the parti- 
tion function (generating function) must be scale invari- 
ant under the renormalization group, we easily get the 
following differential equation 



dg 



(A2) 



with Up a shorthand notation for the dimensionless cou- 
plings, {up} = {k, a, . . . } that occur in 7i, and /J^ the 
corresponding /3-f unctions. 



din 



dx 



(A3) 



with X ~ — ln(L)) the scaling variable. In the vicinity of 
the two-channel Kondo fixed point the /^-functions just 
assume their fixed point value, which are just the renor- 
malization group eigenvalues, y^i = d — x^^, with the di- 
mension d = 1, since all operators are local and live in 
time only. Since, for (f)s we have yt = 1/2, in the close 
vicinity of the two-channel Kondo fixed point we obtain 



AD 

One can also easily show that 



^dg dg 

dD dui 



(A4) 



(A5) 



These relations imply that, g{Lu,T, D) is scale invariant, 
and is only a function of uj/D and T/D. Clearly, similar 
equations hold for the correlation functions of all opera- 
tors with dimension 1/2. Furthermore, the above scaling 
property can easily be modified for operators having di- 
mensions 7^ 1/2. 
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